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Abstract 

The critical properties of renormalizable O(N) field models are determined 
by means of the high order (> 18) behaviour of convergent linked cluster series 
on finite temperature lattices. It is shown that those models become weakly 
coupled at the phase transition. The critical exponents agree to those of the 
corresponding superrenormalizable 3-dimensional models. Concerning critical 
amplitudes and subcritical behaviour, corrections induced by renormalizable 
couplings are measurable. 
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Quantum field theories at finite temperature are the target of recent intensive 
studies. This is not at least because they are supposed to allow for a quantitative 
investigation of the physics of the early universe, including early phase transitions, 
as well as of heavy ion collisions. There is increasing effort to develop and apply 
techniques and methods suitable for a serious understanding both of the theories 
themselves and of their direct physical (measurable) implications. For instance, 
standard weak coupling perturbative expansions suffer from various problems re¬ 
lated to severe infrared singularities for theories with zero mass excitations. Various 
proposals exist to resolve or to circumvent those problems. Among those are hard 
thermal loop resummations [(Lj, finite volume scaled perturbation theory [^, |3j] , and 
the techniques of scale dependent flow equations [|4j. || and block spin renormaliza¬ 
tion group [j6]. 

A very promising way is to start with well defined lattice held theories. This 
in turn implies computer aided studies both numerically and analytically. Con¬ 
vergent expansion methods allow for arbitrary precise measurements within their 
convergence domain, limited only by available skill and resources. One of those 
convergent techniques is the (vertex renormalized) linked cluster expansion (LCE) 
(e. g. Bii). It is a series expansion of correlation functions with respect to the 
strength fields at different lattice sites are coupled. The weaker those couplings, the 
higher the disorder of the system. Those techniques have recently been generalized 
to apply properly to held theories on finite temperature lattices 0- 

In this letter we present results obtained by applying the LCE to the high tem¬ 
perature phase of (low energy effective) O(N) symmetric scalar models on the lattice 
with renormalizable quartic interaction. The symmetry restoration for the 0(4) case 
has been investigated earlier in [|TTJ . Those models provide one of the fundamental 


building blocks of clektroweak standard models and serve as effective theories for 
the QCD chiral phase transition. The expansion parameter is the isotropic nearest 
neighbour coupling. The coefficients of susceptibility series depend on the quartic 
coupling in a non-trival way. They are uniform in sign by inspection, which implies 
convergence up to the critical temperature. Decoding the high order behaviour, 
detailed information on critical properties is obtained, such as numbers or bounds 
on universal critical exponents and amplitudes. 

The standard canonical form of the action reads 


S(t. k. A) = £ (t 2 + A(4> 2 - l) 2 ) - 2 k Y. ' *(»). (!) 

x£A <xy> 

where the last sum is over nearest neighbour lattice sites on the hypercube A = 
Lq x oo 3 with periodic boundary conditions imposed, and 0 < A < oo. We are 
interested in physical (renormalized) coupling constants as defined by the large scale 
behaviour of the vertex correlations 

Fab{p, ~p) = -^-(m 2 R + p 2 + 0{p 4 )) 5 a , b asp = (p 0 = 0,p -»■ 0), (2) 

x/i 
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and 


rj£rf(pi = • • • = P 4 = o) = fc 4 («, b , c, d), (3) 

fij-agbi = ■ • ■ = Pe = 0) = ^^(ai, ■ ■ ■, a 6 ). (4) 

The C n denote the 0(N)-invariant tensors, C^a, 6) = 5 a ,b, and 

2 n 

C'2n(0 J l'> • • • 5 ^2n) ^ ^ ^ai,ai ^ 271 - 2 (^ 2 , ■■■, a*, ■■■, ®2n), n > 2, (5) 

i=2 


where (~) implies omission. In terms of the standard (connected) susceptibilities, 

C'2n\P , li • • ■ , Cl2n) 


(2n- 1)!! 


X2n = 51 < ^ ai (0)^ a2 (x2)---^a 2 A X 2n) >° 


X2,...,X2n 


^ = Y, < $i(°)$i(x) > c 


( 6 ) 

(7) 


x \i= 1 


(spatial weight), they are given by 



-V (xo-ioh' 

1*2 J \ X2, 


( 8 ) 

(9) 

( 10 ) 


It is the set of correlations such as ©,(0) that LCEs are applied to in the first 
place. We do not give technical or graph theoretical details here. For those the 
interested reader is referred to |L(J- One major point to be noticed that goes far 
beyond a purely graph theoretical rearrangement is the introduction of ” 1-particle 
irreducible” susceptibilities 1 ]. It is not just that the series representation of the latter 
are constructed explicitly. Their explicit use allows for a more precise measurement 
of renormalized coupling constants. This is particularly useful for the computation of 
couplings that are related to \n with n > 2, such as are gn and h R in the present case. 
They normally suffer from the problem that their critical exponents are differences 
of those of the higher Xn of about the same magnitude. This adds to the equally 
severe problem that at fixed order increasing n implies decreasing accuracy due to 
potentiated singularities. Fortunately, having the (LCE) 1PI susceptibility series at 
hand, we can avoid this problem by defining appropriate ratios and expressing the 
latter in terms of the 1PI series directly. 

1 They are to be distinguished from the vertex functions leading to (2j)-([|), cf. eg. &0- 
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To this end, we first notice that connected and 1PI susceptibilities are related 


X 2 = 


H 2 — 


1 — 16kx2 F 


lPI,mod 

!-h 


(1 — 16nxl P1 ) 2 ' 


X4 = 


X6 = 


(l-ie n, x rr 

i 

(1 - 16kx2 PI ) 6 


lPI.mod 
X6 i 


where 


lPI.mod 

Xe 


= ^ PI + 12„(x FPI ) 2 , 

= xr+i2o K -hi 


^ 1 - 16/^ PI 

In turn, we dehne the following intermediate functions. 

red _ J_ X2 _ 1 ^ , re d (_ 

Xl ~<k x?' i-i6K X l n - 2 V 

red h2 d 2 _ , red / K\~ 2u 


_ ~ /(red f 1 _ 

1PI — ^2 \ 1 


K\~^ 


^If ~ A™ 
V .1PI — “W 

X2 


/ lPI.mod \o j 

(h2 ) 2 “2 

lPI,mod 

X6 


~ -4^ 1 - 


-(7+^s) 


-(7+^h) 


A 6 - , 1PI mod \ 3 7 6 l • 

(/i 2 )% V K c' 

We have indicated the expected leading singular behaviour. Possible multiplicative 
logarithmic corrections are included in the amplitudes Al rcd . From this we infer that 


6xT 

/4 e % 


1 - — 

K^c 


g R = -36xl c&d 2 


h R = -216*^%- 1-- • (24) 

V K c J 

Introducing the lattice spacing a(n, A) such that T = l/(L 0 a) and T —> T c as 
a(/c, A) —> a(/c c , A) for hxed L 0) we get in terms of dimensionfull (physical) quantities 

^- ~ Aa{T-T c )~ Ua ~ v , (25) 

'^ph 

V r2 - A & (T-T c )-“ h . (26) 



The results presented below are obtained from the 2- and 4-point suscecptibility 
series to 18th order (and the 6-point correlation to 16th order) available for the 
complete range of quartic couplings. We expect the accuracy of an nth order com¬ 
putations on the L 0 x oo 3 lattice to be at least comparable to that one of the order 
n — L 0 on the symmetric oo 4 lattice, if in both cases the same methods are used. So 
at least lattices with L 0 = 4, 6 are covered by our computation. 

The determination of the critical couplings has been done along the following 
lines (cf. also Jl2], [L3], [14]). The above functions are obtained as series representations 


in the hopping parameter 2k, with coefficients depending nontrivially on the bare 
coupling constants A. Let us denote by /l the truncation of the series representation 
of 



/( \k)= E c„(A)(2k)“ 

rrif<u 

(27) 

at order L > 0, 

/l(A, k) — E c„(A)(2 K) v . 

m f<v<L 

(28) 

The susceptibilitiy series 
tion, the coefficients are 
closest to the origin lies 
behaviour of /, assumed 

are convergent within a circle of radius 2 k c (A). By inspec- 
of uniform sign, which implies the (strongest) singularity 
on the positive real axis. Using this, the leading critical 
to be given by 


/ ft \ — u 

f ~ A ^1- J as k —a k c —, 

(29) 

is determined as follows. The critical coupling k c (A) where the 2-point susceptibility 
diverges is computed as the smallest positive root of o^, using that XjE'E) = 
lim K ^ Kc (A)_ y 4PI (A, k) stays finite at k c . That is, 


o' 1 

A 

O 

se 

(30) 


k c ( A, L ) — k c (A) + — + o(L 4 ). 

(31) 


For the regression, we always choose L > 9. The values obtained agree within the 
error bars to those ones got by the standard ratio criterions on the series representa¬ 
tion of X‘ 2 - In this way, however, the error bars themselves are considerably smaller. 
A possible antiferromagnetic singularity at —k c (A) is shifted to — oo by changing to 
the new variables z = , 2 U , leading to 

K c( A) 

Jl(\z) = '52’CvWz v - (32) 

v<L 

Knowing the critical point k c or z c , the critical exponent u is obtained by the large 
order v behaviour of the ratio 

r„ = = - (i + ^^ + rJ) , (33) 

< 4,-1 " / 
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where R v = o(u 4 ) as v —> oo. The decay strength of R v is sensitive to the presence 


of multiplicative logarithmic corrections. A constant amplitude A as above, (29), 
implies that the large order behaviour of R u is determined by the subleading power 
like behaviour on the boundary of the convergence domain, R u = 0(z/ -1 ~ A ), where 
typically 1/2 < A < 1 [|13|, |15[. On the other hand, if A diverges or approaches 0 
logarithmically as k approaches k c , 


A 


log |1 - — 


(34) 


with r / 0, the decay is much weaker, R u = = 0(^A_). Direct measure¬ 

ment of such logarithmic exponents can hardly be done from the LCE series. But 
the y 2 /df can strongly hints this property. 

The critical amplitude is obtained by comparing the series of / to the right hand 
side of (f29[) and linear regression A u = A + 0{y~ l ). 

We don’t use Pade’ approximants. They give hints to the values of critical 
exponents, but the error bars never happened to be smaller. 

Tab. [l] and show the resulting critical line and the critical exponent 7 both 
for the 1-component and the 0(4) model, as obtained on the 4 x oo 3 lattice. The 
coupling A as defined by A = — ^ or O(N), with 

/cF^e-^+A^ 2 - 1 ) 2 ) 

l ’ 2n ~ ’ 35 


is monotonically increasing with A and 0 < A < 1 as 0 < A < oo. Beyond the 
gaussian (zero temperature) hxed point at A = 0 there is def ini te non-mean held 
behaviour at larger couplings. We identify the finite temperature phase transition as 
the strong coupling part, including a range A min < A < 1 where the critical exponents 
are stable. Tab. |3] arid [4] summarize the critical data on the exponents and amplitudes 
for N = 1,..., 4, obtained from the high T part of the phase diagram. They 
are compared to the data of the corresponding superrenormalizable 3-dimensional 
models with the numbers that we have determined along the same lines as outlined 
above|. 

The results can be summarized as follows. The renormalized quartic coupling 
vanishes at T c proportional to the mass (u g = —v)- Within bounds, the 6 -point 
coupling has a vanishing critical exponent (coh = 0). The most precise data for 
the critical exponents 7 and u are obtained by (|33|) with R v = 0 (—2—), indicating 
multiplicative logarithmic corrections to (^9[). All the critical exponents, i.e. 7 , v as 
well as u g , Uh agree in value to those of the 3-dimensional O(N) <h 4 models. The 
numbers on the critical couplings and exponents are the most precise ones as yet 

2 Cf. also e.g. |jl6[ [Tt] . For an extended list of references on 3-dimensional critical numbers both 

for O(N) field models and N-vector spin models cf. || and Q. For the latter, the best of those 
data are given in |l£| . 
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obtained. They are in reasonable agreement to earlier investigations both by Monte 
Carlo methods [|ll|] and scale dependent flow equations 0. 

For a ) g = — v and c Oh = 0, the ansatz (|29|) yields values on the critical ampli¬ 


tudes (|25|) , (p6| ) summarized in Tab. [4j. Although the error bars are rather large on 
those 3rd generation quantities (after the critical line and the critical exponents), 
we measure deviations of the 3d data. This is in conformity with the observation 
that multiplicative logarithmic corrections are present that are not provided by the 
3d superrenormalizable models, ffence we have a twofold strong evidence that a 
quantitative description beyond the leading critical behaviour by a 3d effective the¬ 
ory needs non-superrenormalizable couplings such as <f> 6 , regardless that the quartic 
coupling vanishes at T c . This might become important outside of the critical region, 
and in particular for first order transitions. 

All the 4d data presented above are obtained on the 4 x oo 3 lattice. Consistent 
results are obtained on the larger 6 x oo 3 hypercube. Yet the relative error bars 
become about twice as large due to a reduction of the effective order by two. As yet, 
the data for temporal extent Lq > 8 are not very predictive to 18th order of com¬ 
putations. There are future ways to cure this. First, the construction algorithms 
are that efficient that the order of the susceptibility series can still be increased 
within available resources. Second, zero temperature ’’background” can be sub¬ 
tracted. Third, convergence behaviour can be improved by including ’’irrelevant” 
next to nearest neighbour interactions. 
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Tabic 1: The critical line k c (A) for the 1-component <f> 4 model on the 4 x oo 3 lattice. 
Also shown is the critical exponent 7(A), that exhibits the high temperature phase 
transition as the strong coupling range of the phase diagram. 


A 

A 


7 

1 

00 

0.075773(21) 

1.2407(83) 

0.95 

6.92952 

0.085272(24) 

1.2392(87) 

0.9 

4.3303 

0.093339(27) 

1.2406(102) 

0.85 

3.25072 

0.10085(4) 

1.2395(106) 

0.8 

2.5836 

0.10789(4) 

1.2372(97) 

0.7 

1.7320 

0.12053(4) 

1.2387(110) 

0.6 

1.1769 

0.13094(4) 

1.1977(85) 

0.5 

0.77841 

0.13845(3) 

1.1731(67) 

0.4 

0.48548 

0.14244(3) 

1.1570(58) 

0.3 

0.27538 

0.14239(3) 

1.1323(65) 

0.2 

0.13418 

0.13860(2) 

1.1007(44) 

0.1 

0.04770 

0.13225(2) 

1.0538(46) 

0.05 

0.01997 

0.12866(1) 

1.0311(24) 

0.01 

0.003457 

0.12573(1) 

1.0066(24) 


Table 2: The same as Tab. 1, but for the 0(4) symmetric model. 


A 

A 

HZ C 

7 

1 

00 

0.30967(16) 

1.4473(202) 

0.975 

39.4673 

0.30800(18) 

1.4516(229) 

0.95 

19.4314 

0.30614(19) 

1.4417(243) 

0.925 

12.725 

0.30418(20) 

1.4310(258) 

0.9 

9.34653 

0.30211(18) 

1.4211(233) 

0.8 

4.12700 

0.29204(17) 

1.3780(236) 

0.7 

2.24446 

0.27772(15) 

1.3380(206) 

0.6 

1.25586 

0.25802(13) 

1.2982(197) 

0.5 

0.67905 

0.23313(11) 

1.2569(184) 

0.4 

0.34514 

0.20532(8) 

1.2164(153) 

0.3 

0.16348 

0.17848(6) 

1.1728(123) 

0.2 

0.07041 

0.15583(4) 

1.1247(101) 

0.1 

0.02370 

0.13826(2) 

1.0696(56) 

0.05 

0.00989 

0.13116(1) 

1.0346(25) 
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Table 3: The critical exponents of the <f> 4 O(N) models for N = 1,..., 4, both in 4 
dimensions (upper values) and in 3 dimensions at zero temperature (lower values). 


N 

^min 

7 

V 

vr\ 

UJg 

Uh 

1 

0.85 

0.90 

1.2400(87) 

1.2406(36) 

0.6300(49) 

0.6301(18) 

0.0193(132) 

0.0183(52) 

-0.6296(95) 

-0.6300(42) 

-0.0016(136) 

-0.0008(72) 

2 

0.90 

0.90 

1.3238(139) 

1.3250(52) 

0.6694(64) 

0.6734(28) 

0.0206(150) 

0.0199(76) 

-0.6683(131) 

-0.6788(75) 

0.0030(191) 

0.0040(126) 

3 

0.90 

0.90 

1.4032(156) 

1.4029(85) 

0.7167(76) 

0.7131(40) 

0.0246(169) 

0.0232(87) 

-0.7132(164) 

-0.7165(115) 

0.0010(262) 

-0.0167(218) 

4 

0.95 

0.90 

1.4469(225) 

1.4504(113) 

0.7356(93) 

0.7361(68) 

0.0257(200) 

0.0213(123) 

-0.7285(276) 

-0.7268(134) 

0.0075(397) 

-0.0104(247) 


Table 4: The corresponding critical amplitudes A 4 and for ui g = — v and 
u>h = 0, both in 4 dimensions (upper values) and in 3 dimensions (lower values). 
The numbers hold under the assumption that no multiplicative logarithmic factors 
occur. 


N 

*4.4 

Aq 

1 

16.61 ± 1.46 
23.72 ± 1.49 

590 ± 103 

1386 ± 174 

2 

15.39 ± 2.24 

20.33 ± 1.76 

464 ± 121 

913 ± 177 

3 

14.99 ± 2.40 

17.87 ± 1.68 

406 ± 114 
672 ± 142 

4 

13.83 ± 2.37 

15.29 ± 1.56 

293 ± 94 
447 ± 105 
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